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ABSTRACT 

Context. Lots of information on solar-like oscillations in red giants has been obtained thanks to observations with CoRoT and Kepler 
space telescopes. Data on dipolar modes appear most interesting. 

Aims. We study properties of dipolar oscillations in luminous red giant to explain mechanism of mode trapping in the convective 
envelope and to assess what may be learned from the new data. 

Methods. Equations for adiabatic oscillations are solved by numerical integration down to the bottom of convective envelope, where 
the boundary condition is applied. The condition is based on asymptotic decomposition of the fourth order system into components 
describing a running wave and a nearly unform shift of radiative core. 

Results. If the luminosity of a red giant is sufficiently high (L > lOOL© at M = 2Mq), the dipolar modes become effectively trapped in 
the acoustic cavity, which covers the outer part of convective envelope. Energy loss caused by gravity wave emission at the envelope 
base is a secondary or negligible source of damping. Frequencies are insensitive to structure of the deep interior. 
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1. Introduction 

Observations with the space-borne telescopes CoRoT and 
Kepler brought wealth of information about solar-like oscilla- 
tions in stars along the red giants branch (RGB) (Mosser et al. 
201 1, Hekker et al., 201 1, and references therein). In his recent 

I review, Christensen-Dalsgaard (201 1) describes these objects as 
"the most interesting areas of space asteroseismology". Here we 
focus on dipolar modes in these stars not only because such 
modes are most important for stellar, in particular red giants (see 
Montalban et al. 201 1), diagnostic but also because they present 
some interesting problems for stellar pulsation theory. 

What makes dipolar oscillations different from those of 

• higher degrees is that the Cowling approximation is not appli- 
cable. This approximation has been the basis of standard mode 
classification and analysis of wave propagation. The problem of 

' a unique dipolar mode classification was solved by Takata (2005, 
2006b). The aspect of our interest in this paper is the effect of the 
wave trapping in the acoustic cavity which is essential for under- 
standing oscillations observed in red giants. 

The gross picture of red giant oscillation emerging from 
analysis of massive data on red giant oscillations (Mosser et al. 
2010) turns out rather simple. The power spectra look similar to 
that of the Sun from spatially unresolved measurements. The se- 
quences of low degree modes, up to = 3, are very clearly seen. 
Only after a closer look a difference in dipolar mode pattern may 
be noticed. This difference is important and interesting. 

Observations largely confirmed what has been predicted by 
theoretical studies of Dupret et al. (2009) and Montalban et al. 
(2010). The first paper was devoted to stochastic excitation of 
oscillation in red giants and gives values of expected ampli- 
tudes and lifetimes for five selected red giants models. The sec- 
ond paper presented an extensive survey of the trapped mode 
frequencies in red giant models. In very dense eigenfrequency 



spectra, the modes were identified as the ones corresponding 
to local minima of inertia. Such choice of the modes was mo- 
tivated by the results of the first paper and two earlier papers 
(Dziembowski et al, 2001, Christensen-Dalsgaard, 2004), where 
importance of mode trapping phenomenon in red giant oscilla- 
tions was stressed. The goal of the present work is to improve 
understanding of this phenomenon and its consequences for red 
giant seismology. 

The next section is devoted to derivation of an asymptotic 
approximation for dipolar oscillations in the gravity propaga- 
tion zone and a discussion of transition from standing to running 
waves limit. Based on that, in Sect. 3, the inner boundary condi- 
tion for the acoustic oscillation of the convective envelope is for- 
mulated. Properties of the acoustic cavity in red giant envelopes 
are reviewed in Sect. 4. Efficiency of trapping radial and dipolar 
modes is discussed in Sect. 5. Numerical results on mode fre- 
quencies and damping rates caused by wave losses in selected 
envelope models are presented, respectively, in Sect. 6 and 7. 
The final section confronts theoretical results obtained in this 
and earlier papers with observations. 



2. An asymptotic approximation for dipolar modes 

We consider adiabatic oscillation of an isolated star ignoring all 
effects of rotation. The linearized momentum equation, which 
we write in the standard notation adopted e.g. by Aerts et al. 
(2010) 

(Jp^ = Vp' + pV(D' + p'gCr, (1) 

combined with the Poisson and continuity equations leads to the 
familiar fourth order differential equations for the radial eigen- 
functions. 
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That the system reduces to the second order at{-0 (radial 
modes) has been known from the nascent of the theory of stellar 
oscillations but that this is true also at { = I (dipolar modes) 
has been demonstrated not long ago by Takata (2005). Here we 
we will use his formalism and mostly his notation. The four ba- 
sic radial eigenfunctions ykir) used in this work are defined in 
expressions for vectorial displacement 

^ = r\yrir)er + yiirWH^T^-''', (2) 
and the Eulerian perturbation of the gravitational potential 
VO)' = g\y,ir)er + y3(0VH]Ff e"'"'. (3) 
Note that we have 

where 



B = 



2 

CO r 



We will use standard dimensionless coefficients calculated from 
stellar models 



V =^ 



dlnp 4npr^ 
A = — - Vg, and U = 



dlnr 



Mr 



Takata showed that the radial eigenfunctions of dipolar modes 
satisfy the linear algebraic relation 



UB(y2 -yi) + (2-U + B)b + (1 - B)>'4 = 



(4) 



and used it to reduce the eigenvalue problem to two first order 
equations for radial eigenfunctions 



zi =yi + 



and 



Z2=y2 + 



1 -B 



(5) 



(6) 



l-B 

Here, we write these equations in the following form 
drUij \l-A/B + S A-2-S jUai' ^' 



where 



S = 



UB 



(1-5)2- 



For the radial eigenfunctions describing perturbation we 
have from the definition in Eq. 3 



= (I -U)y3+ y4 

dr 

and from Eq. 4 

iU-2-B)y3 + UBizi-Z2) 



JA = 



1 -B 



(8) 



(9) 



At the base of red giant envelope, B oc r" jM^ may be very small 
even for high frequency modes. Let us note that Eq. 7 in the B 
limit has the same form as in the Cowling approximation but 
with modified definition of variables. Nearly the whole radiative 
core, where A > 0, is a region of gravity wave propagation. Once 
the wave number, -s/lA/Br^, is sufficientiy large an asymptotic 
approximation is applicable. 



For any ^ > in the limit AjB :s> 1 , the fourth order sys- 
tem may be asymptotically decomposed into two second order 
equations yielding one solution that varies rapidly with r and de- 
scribes a gravity wave and one that varies slowly and describes 
a quasi-static response to an external perturbation of the gravita- 
tional potential (Dziembowski, 1971). In 1971, the goal of such 
a decomposition was not only to help understanding physics of 
nonradial oscillations in evolved stars but also to use in the inte- 
rior the analytic form of the rapidly varying solution which was 
difficult to calculate numerically. 

In that work no special attention was paid to dipolar modes 
for which also the slow solution may be derived in an analyti- 
cal form. With Takata's equations, the solution for the potential 
which is divergent at r = is eliminated automatically. A su- 
perposition of the slow and rapid solutions will be used in the 
next section to formulate inner boundary conditions for dipolar 
modes trapped in the convective envelope. 

The slowly varying solution of Eqs. 7 to 9, which is easily 
obtained by setting zi= Z2= 0, is 



V ys,A ) 



c. 



1 

1 

-l+B 
2-U + B 



(10) 



where Cj is a constant. This solution is exact and it is equivalent 
to one given by Takata (2006a, Eq. 116). It describes uniform 
shift of star interior in response to the external potential change, 
which in our application is caused by oscillations in the enve- 
lope. 

To obtain the rapidly varying solution, we adopt 



dzk 
dr 



\krZk\, 



where 




(11) 

is the radial component of gravity wave number. Setting in Eq. 7 



we obtain 



dr2 



(12) 



(13) 



(14) 



The asymptotic solution of this equation is 
Zi=Cr Kdr, 

ykr Jro 

where C,- is another constant. If there is no convective core ro = 
and this form is obtained by fitting to the spherical Bessel func- 
tion ji{r/R), which is the solution of Eq. 7 valid near the center. 
The same form is valid for stars with convective cores but then 
To = Tc, which is the core radius. The first and second exponen- 
tial components describe, respectively, outwardly and inwardly 
running gravity waves. At B <s; 1, Eqs. 7 and 14 imply 



Z2 



y,)zi 



2r3^' 



where 

_ dZi A ■ 
Z2 = r^- + 



dr 



-Zi. 



(15) 



(16) 
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With the same accuracy, Eqs. 5 and 6 imply yi = zi and - Zi- 
Then, Eqs. 8 and 9 yield - -UBy\l2 and y^ = {U - 2)yi + 
UB{z\ - zi)- Thus, for the rapidly varying solution, we get the 
following asymptotic approximation. 



( yr,i ^ 

yr,2 
yr,3 



2Zi 

-UBZi 
-UB[Z2 + iU-4)Zi] 



(17) 



If energy losses in the interior are negligible, the amplitudes 
of the inwardly and outwardly running waves are equal and the 
solution is a standing wave. In luminous red giants even at fre- 
quencies as high as the acoustic cut-off the phase at the top grav- 
ity propagation zone is large. Thus, all nonradial modes have 
mixed character The distance between consecutive modes of the 
same ^ > 1 is much less than the distance between consecu- 
tive radial modes (the large frequency separation) and decreases 
during evolution along the red giant branch. The relative ampli- 
tude in the interior is modulated in the rhythm determined by 
the acoustic frequencies of the envelope. With increasing lumi- 
nosity, the amphtude of modulation increases which leads to the 
appearance of modes trapped in the envelope. This happens first 
for modes of higher degrees but ultimately also for { = I. Such 
modes are of our interest. 

The most important evolutionary efi'ect on nonadiabatic 
properties of nonradial oscillations is the onset of significant 
damping in the inner part of the radiative interior (Dziembowski, 
1971). In the first works on excitation of nonradial modes in gi- 
ants (Dziembowski 1977, Osaki 1977), only unfitted envelope 
models were considered. Solutions of the nonadiabatic problem 
were obtained upon neglecting perturbation of the gravitational 
potential (the Cowling approximation). One needed boundary 
condition was obtained by neglecting the reflected wave. The 
justification was that radiative damping in the deep interior 
causes that the reflected wave returns with much reduced ampli- 
tude and may be ignored in some cases, such as low order modes 
of moderate degrees in Cepheids and RR Lyrae stars considered 
in those papers. The modes were found unstable in spite of large 
damping in the interior because the amphtudes there were very 
low as a result of efficient trapping in the envelope. For luminous 
red giants, such a situation was found also in the case of dipolar 
modes (Dziembowski et al. 2001). 

Complex eigenfrequencies determined with the running 
wave boundary condition set in the envelope yield approximate 
values for true stellar eigenfrequencies if |1'/(rft)| » 1, where 
subscript / denotes the imaginary part of a complex quantity (we 
will use subscript R to denote the real part) not only for unsta- 
ble modes but also for stochastically driven stable modes. The 
requirement is that the damping rate is not too high. The explicit 
condition is easy to write in the case when the quasi-adiabatic 
approximation is still applicable and |w/| «; wr. In such a spe- 
cial but not unrealistic case we have (e.g. Dziembowski et al. 
2001) 

kr^ [I -1(0+ )], 



(Or 



where 



D = 



0}R 



(Vad - V), 



and this condition is 



L 



where *Pi - *Ps(rft). The trapping effect follows essentially from 
adiabatic properties of stellar oscillation but it is enhanced by 
nonadiabatic effects in the case of unstable or weakly stable 
modes. 

As long as we stay within the domain of linear oscillation 
theory, the use of running wave solution in the situation when 
the strong inequality in Eq. 18 is not satisfied means an over- 
estimate of damping in the core. However, it must be kept in 
mind that nonlinear effects may significantly enhance radiative 
damping. Could such effects be of any importance in the case of 
stochastically driven oscillations? Dupret et al. (2009) found that 
everywhere in their models |(f • V)^| « |f | and gave the negative 
answer. However, their argument refers to a strong nonlinear- 
ity, whereas at much lower amplitudes the resonant excitation 
of high degree modes may occur which leads to a significant 
increase of the damping rates (e.g. Kumar & Goodman, 1996). 
Unfortunately, in this case, a precise determination of the condi- 
tion for the running wave approximation is not easy. 

In this paper, we adopt the running wave boundary condition 
at the top of radiative interior. This means that we treat the con- 
vective envelope as a separate oscillating body loosing part of its 
energy by emitting gravity wave from its bottom. We will exam- 
ine implication of this assumption and appeal to observations to 
verify its applicability. 

3. The boundary conditions and eigenfrequencies 

The asymptotic approximation developed in the previous section 
allows to formulate the inner boundary condition which may 
be imposed on numerical solutions of the equations for dipo- 
lar oscillation in the outer layers. Radiative dissipation arises 
mainly deep in the star interior, so that even when it is large, 
we may still rely on adiabatic approximation in the outer part of 
the g-mode propagation zone, where we want to apply the inner 
boundary condition. This is not an essential approximation and 
in fact it was not used in Dziembowski (1977) paper. However, 
we will avoid the Cowling approximation, which may be glob- 
ally justified at higher degrees but not, as Christensen-Dalsgaard 
& Gough (2001) stressed, at ^ = 1. We are now going to correct 
this mistake by taking into account perturbation of the gravita- 
tional potential and the slowly varying solution. 

The general solution of Eqs. 7 to 9 is a linear superposition 
of the slow and rapid solutions given in Eqs. 10 and 17, respec- 
tively. It involves two unknown constants and thus yields two in- 
dependent boundary conditions. The slow solution has the form 
specific for dipolar modes and only for these modes it is inde- 
pendent of the interior structure. The form of the rapid solution 
is valid for any { if in Eq. 1 1 for kr the number 2 is replaced with 
{{{ H- 1). For the running wave solution, we get from Eqs. 14 and 
16 



Zi 



irkr + ^\A- 



r d|k,| 
\kr\ dr 



Zi 



(19) 



d^R » 1, 



(18) 



The boundary condition obtained in this way may be applied 
at any place within the radiative interior, where the asymptotic 
approximation is valid. The explicit form of the inner bound- 
ary condition for the case of running wave is given later in this 
section. The same conditions may be applied at the base of con- 
vective envelope, r - r^, if B(rh) « 1 and there is a jump in A to 
a finite value due to discontinuity in hydrogen abundance. Note 
that then a corresponding discontinuity of y2 must be taken into 
account. However, it is more realistic to treat A{r) as a steep but 
continuous function. In such a case, some modifications in the 
form of boundary are needed. 
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To this end, we consider thin layer beneath the envelope 
bottom where A rises from zero at r = r/, to a finite value at 
its lower boundary where A/B becomes large. Like before, we 
start with Eq. 13 but now we seek its asymptotic solution in the 
B — > limit in terms of the Airy function and Si (see e.g. 
Eq. 10.4.112 in Abramovitz& Stegun 1972). The solution is 



where 

s = -A^'\n - r), A 



(20) 



2 dA 
r^B d7 



> 0. 



Note that 



dr 



dZi 
ds 



where 



2rdA\ 
B d7j. 



1/3 



(21) 



At (rh - r)K » rb, 
Zi oc|ir'^'*(C4sin>P + CBCOs»F), 



3^-^^^ ^ 4- 



Thus, for the inwardly running wave we must have Ca — -iCr, 
Cr = Cb, and we get 

Zi = Siis) - ims). (22) 

The expression for Z2 in Eq. 15 remains valid but with 

d A-V„ 
Zi = /3-iSi - im + — ^Zi (23) 
ds 2 

Again with Eqs. 8 and 9, we may derive the corresponding ex- 
pressions for y3 and ^4. The slow solution given by Eq. (10) re- 
mains unchanged. For the explicit form of the solution at con- 
vective envelope bottom, we use the values of the Airy functions 
and their derivatives at i = 0, 



Vf' 



"ds" 



-1 dSi 
V3 ds ' 



1 dSi 
S/"ds" 



r(2/3) 
r(i/3) 



0.729, 



and A{rh) - 0. 

At r = rh, with the use of Eqs. 10, 17, 22, and 23, we finally 
obtain 



yi 
y^ 

yy4J 



1 \ 
1 

-l+B 
2-U + B 



2Zl 

Zi 
-UBZi 
-UB[Z2 + (U-4-)Zi] 



(24) 



where 



Zi = 1 



rh 



Cr 



1 

Vf' 



and 



Z2 = 0.729/? 1 



V3 



'-Zi. 



The same form is valid for the inner boundary condition applied 
within the radiative core but, as follows from Eq. 19, with Zi - 
1, and 



Zi = UK + 



Zi 



Two inner boundary conditions in the form of linear alge- 
braic relations are obtained by eliminating constants C, and Cr- 



The relations are complex and so are the resulting eigenvalues 
and eigenfunctions even if, as we will do in this paper, all nona- 
diabatic effects in the envelope are ignored. 

Derivation of the outer boundary conditions for oscillations 
in realistic stellar models cannot be done in such a rigorous way 
as in the case of polytropic models considered by Takata (2005). 
In realistic models, Vg(R) < 00, so there is no singularity, and 
t/ > 0. In main sequence stars similar to the Sun, V;, ~ A ~ 10^ 
and U ~ 10 In stars near the tip of the RGB, ~ A ~ 10^ 
and U ~ 10"^. The standard procedure is to set the boundary 
condition in the atmosphere near maximum of Vg and ignore U. 
Then in an equation like Eq. 7, the coeflicients are approximately 
constant and for the wave reflected downward we have 



l-Il 

2 



1- A 1- 



yu 



(25) 



where oj^c = g^i/^c is the acoustic frequency. The vacuum so- 
lution of the Poisson equation outside the star implies 3^4 = -2y3 
at r = R. This is the second outer boundary condition adopted 
in this paper For dipolar modes, an additional condition fol- 
lows from Eq. 4, which combined with the previous two yields 
^3 - U(yi - 3'2)/3. Thus, we have three conditions and, with 
the standard normalization yi - 1, we may carry a single inte- 
gration from the surface downward. Application of this standard 
procedure to luminous red giants may seem problematical and 
therefore the sensitivity of results to variation in the outer bound- 
ary condition must be checked. Only one of the inner boundary 
conditions is needed in a search for eigenfrequency. The second 
condition may be used as a test of numerical accuracy. 



4. The acoustic propagation zone 

It is well-known that frequencies of low order radial modes in 
evolved stars may be accurately calculated from their envelope 
models alone. We will see that in the case of luminous red gi- 
ants, this property extends up to the highest orders. The bottom 
of the envelope, which must cover the whole acoustic propa- 
gation zone, may still be within the convective region. Thus, 
what determines the depth of the propagation zone is not the 
N{r) = a> condition, where N is the Brunt-Vaisala frequency. 
A general condition limiting the propagation zone from its bot- 
tom does not exisfl The condition a> < for wave reflection 
is derived for an isothermal plane-parallel atmosphere built of 
perfect gas. Although in this case the condition is numerically 
close to w < TV, its origin is different (to see this consider, e.g., 
the case of gas with Fi = 1 when N = 0). In fact, the condi- 
tion follows from general rule prohibiting wave propagation if 
the background structure changes significantly over wavelength. 
The equality ll) - co^c corresponds to equality of the radial wave 
number acoustic waves to 0.5 /Hp, where Hp = p/gp is the 
pressure distance scale. It does not seem unreasonable to use 
o) < (jJacif) as a crude local no-propagation condition, as long as 
Hp is the shortest of all local distance scales. 

For nonradial modes the acoustic propagation zone may 
be limited from the bottom by the Lamb frequency, £.( - 
+ l)/r exceeding u. This critical frequency was determined 
upon assuming the Cowling approximation. Avoiding this ap- 
proximation, Takata (2006b) found a simple modification of the 



' The generalized acoustic cut-off frequency, u)c, as defined e.g. in 
(Aerts et al. 2010) is not a good choice in the present application be- 
cause (jj^ < in the lower part of the envelope. 



4 



W. A. Dziembowski: Dipolar Modes in Luminous Red Giants 



Table 1. Selected models from the BaSTI track at Mq = 2Mq 
andZ = 0.02. 





M 


logL 


lOgTeff 


n/R 


Vmax 


Av 




(Me) 




(K) 




(pHz) 


OuHz) 


1 


1.998 


1.500 


3.692 


0.0746 


110.8 


8.851 


2 


1.997 


2.000 


3.663 


0.0454 


27.47 


3.071 


3 


1.994 


2.501 


3.628 


0.0287 


6.51 


1.050 


4 


1.985 


3.001 


3.586 


0.0104 


1.46 


0.329 


5 


1.967 


3.355 


3.552 


0.0061 


0.49 


0.137 


6 


1.967 


1.771 


3.687 


0.2410 


55.90 


5.752 



Lamb frequency for dipolar modes, which in our notation may 
be written as 

x>™ = x.(i-f). 

In our appUcation, this modification is not very significant. Near 
the bottom of convective envelope the correction is at most one 
percent. The local maximum of f//3 (up to 0.5) is reached in the 
middle of the envelope. We will not use modified N, where the 
same factor appears in the denominator, because it has no efi'ect 
on lower boundary of the acoustic propagation zone. 

During evolution along the RGB frequencies of stochas- 
tically excited oscillations decrease. The critical frequencies 
<^ac('') and JLimir) at specified fractional radius decrease too, 
but not as fast. Therefore the bottom of the acoustic propaga- 
tion zone moves upward. In the same time, the bottom of the 
convective envelope moves downward. This leads to expansion 
of the evanescent zone between the gravity and acoustic prop- 
agation zones and to efficient mode trapping in the convective 
envelope. 

As an example, we consider a sequence of red giant enve- 
lope models with mass, M, luminosity, L, and efifective temper- 
ature, Teff, taken from BaTSI (Pietrinferni et al. 2006) evolu- 
tionary tracks for the initial mass Mq = 2Mq and heavy ele- 
ment abundance parameter Zq - 0.02. The opacity and equation 
of state data come from the same source. The adopted mixing 
length parameter a = 1 .85 corresponds to 1 .74 in BaTSI mod- 
els. Parameters of these models are listed in Table 1 . The models 
may be regarded realistic down to the bottom of the convective 
envelope. Listed in the table values of Vmax, corresponding to the 
expected maximum of acoustic power for stochastic oscillations, 
were evaluated with the Kjeldsen & Bedding (1995) expression 

M I T \^-^ 

Vmax = Y| — I X3050/Z//Z. (26) 

The values of Av, which denote the distances between consecu- 
tive radial modes (the large separations), were calculated around 
Vmax for reconstructed envelope models with Warsaw codes (see 
e.g. Dziembowski & Soszynski 2010). Models 1 to 5 represent 
stars ascending the RGB. Model 6 in Table 1 represents a heUum 
burning (red clump) object 

We begin with the extreme case, that is, with Model 5, which 
is located at the tip of red giant branch. The three critical fre- 
quencies in this model are shown in the top panel of Fig. 1 . The 
acoustic cavity for radial mode is determined solely by the run 
of a>ac. For dipolar modes the bottom of the acoustic propagation 
zone is determined by 

OJ = max[wac('-),-Cim('")] 



condition and the top of the gravity propagation zone hy oj = N. 
In the middle and bottom panels, we show fractional inertia 

m= f (|);i|2 + 2|);2lV'dr (27) 

for radial and dipolar modes, respectively. For these plots, nor- 
malization / = I(R) = 1 was adopted. With the standard normal- 
ization yiiR) = 1, the value oi I = /(/?) is a measure of mode 
trapping. Let us stress that / is not the total mode inertia but 
only its part contributed by the convective envelope. The plots 
may suggest that for any of the modes shown in Fig. 1, there 
is no contribution to the total inertia (energy) from the interior 
below r » 0.2. In fact, for the dipolar modes the inertia, Ic, aris- 
ing from the slow component, which keeps the center of mass of 
the whole star at rest, is substantial. This contribution to the to- 
tal inertia may be calculated without any knowledge of the core 
structure. Using in Eq. 27 the expressions for yi and j2 from Eq. 
10, we get 

Ic = \l>'ul' + 2|j,.2l V'dr = 3 ICsl^ (28) 

The Ic/I ratio ranges up to 1/3 for the po,i mode and decreases 
with the mode order. 

According to the terminology adopted in works on stellar 
pulsation, the radial pi is called fundamental (no nodes between 
the surface and the lower boundary and one zone of steep rise 
in /). Its frequency is well below the minimum value of cj^dr)- 
For the overtones, most of the contribution to / arises in the zone 
^ (^dcir) condition or close to it. Comparison of middle and 
bottom panels reveals first of all a similarity of the energy dis- 
tribution in dipolar and radial modes, particularly in case of first 
(fundamental) modes whose frequencies are close. More differ- 
ences are seen for the overtones. 

Let us note in the top panel that the critical frequencies coheir) 
and Jlimir) are close to each other in the frequency range of the 
dipolar overtones depicted in the panel. The role of the Lamb 
frequency is reflected in the shape of /(r) for all dipolar over- 
tones. Beginning with p3, this critical frequency determines the 
bottom of the acoustic propagation zone. 

Where between consecutive radial modes the dipolar mode 
appears depends on the critical frequencies determining the 
lower boundary of the acoustic propagation zone. If it is 
which is ^-independent, the dipolar and radial modes have close 
frequencies. This is the case for the low order modes in Model 
5, as we may see in numbers quoted in Fig. 1. With increasing 
oscillation frequency, the role of the Lamb frequency increases 
and the dipolar mode moves away but still is rather far from the 
mid point between consecutive radial modes. It comes closer to 
this point in less luminous stars where higher order modes occur 
near Vmax- 

Fig. 2 shows in separate panels fractional moment of inertia 
for modes with frequencies around Vn,ax in Models 1 to 4. The 
places where the critical frequencies are equal to Vn,ax are marked 
in each panel. The bottom of the acoustic cavity is clearly within 
the convective envelope except for Model 1, where the equality 
Wac = 2;rvmax takcs place in the radiative core. For the second 
critical frequency, the equality occurs in the envelope but the 
evanescent zone is too narrow for an efficient mode trapping of 
the dipolar mode within the convective envelope. We may ob- 
serve in this figure how the width of the evanescent zone in- 
creases with stellar luminosity leading to trapping { = and 1 
modes in the upper part of the envelope. This is seen in the up- 
ward shifts of the oscillatory pattern in I(r). 
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Fig. 1. Critical frequencies (see text) in Model 5 (see Table 1). 
{top) Relative contributions to mode inertia integrated upward 
from the base convective envelope for all four modes with fre- 
quencies below Lo^ciR)- (middle) The same for the first four 
dipolar modes trapped in the envelope, (bottom) There is one 
more dipolar overtone, whose frequency v = 0.677/jHz is below 
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Fig. 2. Relative contributions to mode inertia integrated upward 
from the base convective envelope for three modes (two radial 
and one dipolar) with frequencies close to v,nax in the four in- 
dicated models whose parameters are given in Table 1. The 
short lines parallel to the ordinate axis mark the bottom of 
the convective envelope (long-dashed style), the places where 
2;rymax = (^ac, and Xim (drawn in the solid, and short-dashed 
style, respectively). The styles correspond to those used to show 
the respective frequencies in Fig. 1. 
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Fig. 3. Critical frequencies in Model 6. The two horizontal lines 
correspond to dipolar mode frequencies. The lower one to pi_i 
and the upper to pg 1, whose frequency is near Vn,ax- 



Fig. 3 shows critical frequencies in Model 6 representing an 
object in the phase of He burning in the convective core. The 
convective envelope is relatively shallow and the - Inv^^^ 
equality takes place in the radiative part of the interior. In Sect. 
3, we described the version of the running wave boundary con- 
dition, which may be applied in the gravity wave propagation 
zone. However, it is rather unlikely that the assumption that the 
reflected wave has negligible amplitude might be justified in this 
case. 



5. Mode trapping 

Mode trapping in the acoustic cavity is a result of wave reflection 
within the evanescent zone. The efliciency of this process de- 
pends on the width of this zone and on mode degree. The eigen- 
function in the evanescent zone is a superposition of monotoni- 
cally increasing and decreasing components. At higher degrees, 
the dependence of the displacement on the distance to the cen- 
ter in the two components is approximately given by ' and 
'pjjg jarge contrast between these behaviors causes one of 
the component to dominate and it is the increasing component 
for modes trapped in the envelope. A narrow evanescent zone 
suffices for an efficient trapping. We may speak about hard re- 
flection. This effect is behind the occurrence of unstable modes 
of moderate { in models of Cepheids and RR Lyrae stars. At 
C - Q and 1 there are also monotonically increasing and de- 
creasing components but not as fast as at higher degrees. The 
reflection is softer and a wider evanescent zone is needed for an 
efficient trapping. 

For radial modes, efficient trapping manifests itself in a weak 
dependence of frequencies on the exact form of the inner bound- 
ary condition set at the bottom of the convective envelope. This 
is so because then the contribution from the decreasing compo- 
nent becomes very small once the acoustic cavity is reached and 
thus have a negligible influence on eigenfrequencies. The effect 
may be quantified by considering the generic form of the bound- 
ary condition for homogenous second order equations 



dyi 

cos(7rcb)r h sm(7T(p)yi 

dr 







(29) 



and letting <p vary between and 1. The numbers quoted in Figs. 
1 and 2 were obtained assuming (p - 0. In Fig. 4 the effect is 
illustrated for Models 1 and 2. In the former model, except of a 
narrow range, the value of <p has almost no effect on frequencies. 
In this range, the boundary condition allows a large decreasing 
component which is certainly impossible for any realistic model 
of the whole star. With the yi(ri,) = condition we are safe 
away from this bad range. In Model 1 the effect of varying (p is 
considerably larger which is consistent with what we see in Fig. 
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Fig. 4. Dependence of radial mode frequencies, v, on the param- 
eter (f) adopted in Eq. 29. The pairs of frequencies around Vmax in 
Models 1 and 2 are considered. Plotted are the shifts relative to 
the values at cp - 0. The top panel shows the shift in the lower 
frequency. The bottom panel shows the shift in the frequency 
difference. 



2. In this model, the propagation zone for the two radial modes 
extends to the core. 

Changes in dipolar mode trapping may be observed by relax- 
ing the inner boundary condition and treating /, defined in Eq. 
27, as a continuous function of frequency. The top panel of Fig. 
5 shows the values of / in the frequency range centered at v,nax 
for Models 1 and 2. The inertias of radial modes in the same 
frequency range are shown for a comparison. It is important to 
notice the difference between the two models. The minima for 
Model 2 are far deeper than in Model 1 . The difference in trap- 
ping is better visible in the lower panel showing values of \y i +yT, \ 
at the lower boundary. In this combination of the eigenfunctions, 
the contribution of the slow component in is nearly eliminated, 
as we may see in Eq. 24. In this case the minimum value is zero. 
The dependence is so steep in the case of luminous red giants 
that unfitted envelope models suffice to determine frequencies 
of expected dipolar mode frequencies at specified stellar param- 
eters. 

The transition to efficient trapping occurs at logL « 2. 
Luminosity of our Model 1 places it between Models A and B of 
Dupret et al. (2009) while that of Model 2 between their Models 
B and C. In the latter model (log L - 2.1) they find that damping 
in radiative interior is high enough to reduce amplitudes of all 
dipolar modes below the detection level, except the ones corre- 
sponding to / minima. The value of the critical luminosity refers 
to models with Mq = 2 and it is lower at lower mass. 



6. Frequencies of the trapped modes 

The minima of \yi + y^lt, shown in Fig. 5 occur at frequencies 
close to real part of the eigenfrequencies obtained with the code 
described in Sect. 3. The location of the minima and the dipo- 
lar mode frequencies close to the mid points between consecu- 




Fig. 5. Inertias of dipolar and radial modes in the frequency 
range centered at Vmax in Models 1 and 2 are shown in the 
top panel. For dipolar modes the inner boundary condition was 
relaxed so that inertia is a continuous function of frequency. 
The standard outer boundary condition and normalization were 
adopted in all the cases. The quantity plotted in the bottom panel 
is a measure of the gravity wave amplitude at the top of the prop- 
agation zone. 



tive radial mode frequencies is expected in the B ^ oo limit. 
Note that in this limit 5 — > and Eq. 7 takes again the form 
corresponding to the Cowling approximation, which is the basis 
of the usual asymptotic expressions for dipolar mode frequen- 
cies. However, asymptotic expressions yield a good approxima- 
tion for eigenfrequencies only if at the place where X-iir) = 2nv, 
we have 15 1 <k 1 and this is not true in most the cases considered 
in this paper. Beginning with £ - 2, the Cowling approximation 
is applicable and the asymptotic expressions may be derived in 
the usual way. 

The real parts of eigenvalues obtained in the way described 
in Sect. 3 differ very little from eigenfrequencies obtained with 
the simple real inner boundary condition yKr/,) - 0. Only in 
Model 1 the difference exceeds one percent. 

The small separations based on dipolar and radial mode fre- 
quencies are one of standard characteristics of stochastic oscil- 
lations. Here, we use relative small separations defined as 



SV0,l _ 0.5(Vn+l,0 + V„,o) - Vn,l 

Av 



Vn+1,0 - Vnfi 



(30) 



which measure the departure of dipolar mode frequencies from 
the mid position between consecutive radial modes. In the top 
panels of Figs. 6 and 7 we show values of 6„ for modes at 
y ~ Vmax and n - I, respectively. The values were calculated 
for envelope models with the surface parameters taken from 
BaSTI tracks for metallicity parameter Z - 0.02 and initial 
masses Mq - 1 and 2Mq. In Table 2 we list selected parame- 
ters for selected models at Mo = IMq. Similar data for models 
at Mq = 2Mq and nearly the same L were given earlier in Table 
1. 
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Fig. 6. The normalized small separation defined in Eq. 30 (top) 
and the lifetimes limited only by wave energy loss (bottom) for 
the dipolar mode with frequencies closest to Vmax in red giant 
envelope of indicated luminosity and initial mass. 




Fig. 7. The same as in Fig. 6 but for pi i mode 

Table 2. Selected models from the BaSTI track at Mq - IMq 
and Z = 0.02. 





M 


logL 


log Teff 


rtIR 


^max 


Ay 




(Me) 


(Le) 


(K) 




O^Hz) 


(//Hz) 


1 


0.997 


1.500 


3.658 


0.0628 


41.53 


4.869 


2 


0.994 


2.000 


3.628 


0.0410 


10.31 


1.717 


3 


0.986 


2.501 


3.590 


0.0213 


2.38 


0.586 


4 


0.964 


3.005 


3.542 


0.0125 


0.49 


0.172 


5 


0.933 


3.275 


3.511 


0.0132 


0.20 


0.092 



The jump-Uke changes, visible in the first figure, result from 
gradual switches of mode orders, from n - 8 to 3 and from 
n - 12 to 4 for M - \ and 2Mq, respectively. The changes at 
constant n, as we may see in Fig. 6, are smooth but not nec- 
essarily monotonic. At n = 1 , the dipolar mode occurs always 
significantly closer to the higher frequency radial mode. Clearly, 
negative values of 6„ prevail also for modes near v,nax- Only in 
the low luminosity models at Mq - 2Mq the sign is opposite. 
This agrees with results of Montalban et al. (2010). Negative val- 
ues of 6„ were also found in the fit to the CoRoT data (Mosser 
et al. 2010). Significant departure from the mid position means 
that 6„ is an independent seismic observable, which may per- 
haps be measured with future, more precise data. As for today, 
only mixed dipolar modes give us new seismic information about 
red giants beyond Av and Vmax- We will return to this matter in 
Section 8. 

7. Effect of wave losses on mode lifetime 

The imaginary part of the eigenfrequency calculated with our 
code yields the damping rate, y = -cl>i, which takes into account 
only the effect of the wave losses at the bottom of the convec- 
tive envelope. Of course, this is not the only important source 
of damping but the other resides predominantly in the outermost 
layers, where cj^^ -Ce and the eigenfunctions depend mainly 
on mode frequency and not on its degree. Therefore, the dif- 
ference between dipolar and radial mode damping rates may be 
assessed without explicit calculation of nonadiabatic effects in 
these layers. 

Proceeding a similar way as in derivation of the variational 
expression for real eigenvalues, we begin with multiplying Eq. 
1 by ^' and integrating over the envelope. In further transforma- 
tion, we make use of the Poisson and continuity equations and 
finally arrive at 



n\pf 



■Agp- 



AnG dr 

IVOf 
' AnG 



dl + 



(31) 



where E is the ad hoc added ^-independent term describing the 
near surface effects that may contribute to damping (or driving). 
Assuming <K ojr, we get from the imaginary part of Eq. 

31 



W + E, 
IojrI 



where 



W 



O' (90'* 



AnG dr 



dS, 



which in terms of the radial eigenfunctions defined in Eqs. 2 and 

3 becomes 



W 



^(<^yly2) + -^^(y3yl) 



With the use of the asymptotic solution given in Eq. 24, a more 
revealing expression for the damping rate caused by the wave 
losses may be obtained. Neglecting terms of the orders of B and 
<jDiI<jDr in comparison with 1, we get 



U 
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The leading terms in the coefficient at the C,Cr product cancel 
out. From Eqs. 14 and 16, we have 

3(ZtZ2) = 0.842yS. 

Thus, the total damping rate, y - -coi, is given by 

7 = yw + 7E = -(0i= ^0M2o}R\CA^P(pr\ + ^ j . 

The expression for /<;, which is the part of mode inertia arising 
in the core, is given in Eq. 28. The rapid component contributes 
only to the energy loss because the reflected wave was neglected. 

The near surface effects have some influence on the eigen- 
functions at the bottom of the envelope and, hence, on the value 
of W. This influence was assessed by adding to the coefficent 
at yi in the r.h.s. of Eq. 25 an imaginary number leading to the 
values of ys of the order of those found by Dupret et al. (2009). 
The changes in values of Jh, were found small enough to regard 
the two dominant sources of damping as additive. 

Damping due to wave losses sharply decreases with increas- 
ing red giant luminosity. This is shown in the bottom panels 
of Figs. 6 and 7 for the dipolar modes with frequencies clos- 
est to Vjnax and for pi i modes, respectively. The plots depict, for 
a number of envelope models, the values of tw = 7^, that is, 
lifetimes limited only by the wave emission. There is some un- 
certainty in evaluated lifetimes following from the uncertainty in 
the derivative dA/dr, at the top of radiative interior. This deriva- 
tive enters the inner boundary condition (see Eq. 24) through 
parameter jS defined in Eq. 21. Our standard approach is to eval- 
uate a> and tw assuming constant element abundance but the He 
abundance may rise immediately beneath the bottom of the con- 
vective envelope. To check consequences of this uncertainty, cal- 
culations were performed with the derivative increased by factor 
10. As expected, the effect on frequencies was negUgible. The 
lifetimes became by factor about two longer. This is not much 
for such a drastic modification of the A derivative and, especially, 
when the effect is compared with the four orders of magnitude 
increase of tw between logL = 2 and 3. The general pattern of 
the TwiL) dependence, which we see in these two figures, are 
similar. The lifetimes are longer for the pu modes because at 
lower frequency the evanescent zone is wider, hence trapping is 
more efficient, and the wave amplitudes are lower. Decreasing 
mass at a specified luminosity acts in the same direction. 

Let us now focus on modes with frequencies Vmax which are 
of interest for comparison with observations. The lifetime due 
to wave losses increases rapidly with L and it scales roughly 
as M/L. The average exponent, w, in the tw ~ relation is 
close to 4. It is interesting to compare the role of wave losses 
and the near surface damping in limiting mode lifetimes. For the 
latter effect we use values of te calculated for radial modes of 
by Dupret et al.(2009). They calculated also lifetimes for non- 
radial modes but these depend on mode properties in the core. 
Therefore, it seems safer to rely on results radial modes at simi- 
lar frequencies. For the trapped dipolar modes, which we want to 
consider, Ic/I is small and / is close to mean value of the nearest 
radial modes. So that, we may use mean of the lifetimes of these 
modes as a crude assessment of te for trapped dipolar modes. 

In models with Mq = 2Mq at logL = 1.8 and 2.1, the val- 
ues of Vmax are 49 and ll/Aiz , respectively. From our Fig. 6, 
we find the values 20 and 110 d for tw- The corresponding val- 
ues of Te inferred from Figs. 7 and 8 of Dupret at al. are 20 and 
30 d. Thus, we may conclude that at Mq = 2Mq and logL > 2 
the wave losses become a secondary source of damping and be- 
come negligible at somewhat higher luminosity. At lower ini- 
tial mass the same happens at accordingly lower luminosity. At 



Mo = 2Mq and logL = 2.1 the total relative damping rate, - 

0)i/(Or = (VmaxTT^ = 0.024. 

With this estimate, we may return to the matter of validity of 
the running wave approximation, discussed in Sect. 2. In order to 
use condition given in Eq. 18, we still need the values of *Pi{ and 
J Dd^R for which complete stellar models. For the present esti- 
mate we will use numbers for the M - 2Mq models calculated 
for Dziembowski et al. (2001) paper which should be adequate 
for this aim. The values of "Pr and / Dd^PR at Z, = 2 are 1 .9 x 10^ 
and 0.49, respectively, while at L = 2.5 the corresponding values 
are 1 .0 x lO"' and 370 (in this case the required exact expression 
for kr was used). The value of J Dd*PR grows much more rapidly 
with L than and it is clear that at the value sUghtly exceeding 
2 the running wave approximation is justified even if only the 
linear radiative damping is considered. 



8. Discussion 

Solar-like oscillations in low luminosity red giants are inter- 
preted as global mixed modes with both the gravity and acoustic 
cavities contributing, in various proportions, to the total inertia. 
At higher luminosity, a more adequate picture is that of acous- 
tic modes of the convective envelope loosing energy via gravity 
wave emission from its bottom. With growing luminosity this 
loss becomes insignificant. 

Transition between these two pictures is continuous but very 
sharp. From Dupret et al. (2009) and the arguments presented in 
this paper it follows that for stars of the initial mass Mq = 2Mq 
the transition occurs near luminosity L > 100, which corre- 
sponds (see Table 1) to v„^ax = 27.5 and Av = 3.1/i. At lower 
mass, the transition takes place at lower luminosity but the values 
of Av at the transition should be similar. In Fig. 3 of Mosser et 
al. (2010) most of data points lie above Av = 3/zHz but there are 
also lots of data points below Av = 2//Hz which are most likely 
associated with luminous objects, where only acoustic envelope 
modes are expected. The most striking feature shown in this fig- 
ure is the universal pattern of the relative positions of the dipolar 
and radial modes across the whole range of Av values from 1 
to 8/iHz. The dipolar modes are located between radial modes 
somewhat shifted off the mid point toward high frequency. This 
is consistent with what is seen in Fig. 5 for the envelope modes. 
We may also see that some spread of the position of such modes 
is expected. Similar situation is expected for the mixed dipolar 
modes corresponding to minima of the inertia (Montalban et al. 
2010). 

A convincing evidence for the mixed mode excitation comes 
from detection of few dipolar modes near the mid points be- 
tween consecutive radial modes (Beck et al. 201 la. Bedding, et 
al. 2011, Mosser et al. 2011). In the first two papers, three ob- 
jects from the Kepler catalog are discussed. The values of Av 
are in the range of 7.9 to lO/zHz where stochastic excitation of 
mixed modes is indeed expected. However, Mosser et al. (201 1) 
in their analysis of the CoRoT data find a large number of such 
objects in a wide range of Av extending down to 2 /iHz, that 
is well below the range expected range of mixed modes. They 
also report an evidence for mixed { - 2 mode excitation and 
note this is in conflict with the Dupret el. (2009) prediction that 
such modes should be eUminated by strong damping in the core. 
A possible solution of these two problems could be the gravity 
wave reflection caused by a steep composition gradient above 
the hydrogen burning shell, where most of the damping takes 
place. The large values of AJobs determined by Mosser et al. for 
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the dipolar mixed modes at low Av are consistent with such an Takata, M., 2006b, Pub. Astron. Soc. Japan, 58, 893 
interpretation. 

Acoustic envelope modes are totaly ignorant of what hap- 
pens below the convective envelope. Their frequencies yield in- 
dependent but similar information to radial modes because in 
both cases, frequencies are determined by sound speed distribu- 
tion within the envelope. Only mixed modes probe the radiative 
interiors of red giants. Data on multiple peaks has been already 
used to distinguish between helium burning and other red giants 
(Bedding et al. 2011, Mosser et al. 2011). There is a prospect 
for going beyond that. In particular, we may get interesting con- 
straints on element distribution from frequencies and on dissipa- 
tive processes in the core from amplitudes and lifetimes. A new 
exciting result based on mixed dipolar modes is the evidence for 
fast rotation in red giant cores (Beck et al. 201 lb). 

We do not yet have data satellite on oscillation in most lu- 
minous red giants (logL > 2.5, Av < 1). Rich data on such ob- 
jects located in the Magellanic Clouds are available from OGLE 
project (Soszynski et al. 2009). The advantage of these data is 
that luminosity of individual objects is known to a good accu- 
racy. The peaks occur at frequencies below l/^Hz, which is con- 
sistent with prediction based on Eq. 26. The frequencies are in 
the range of lowest order p-modes. In a number of cases two 
clearly dominant peaks are seen. However, the frequency sepa- 
rations are by some 20% lower than expected large separation. 
This could not be explained in terms of two f = modes, nor in 
terms of the { - 0,1 pair (Dziembowski & Soszynski 2010). In 
the light of the discussion in the present paper the second option 
would appear more plausible but, in fact, it leads to a greater dis- 
crepancy. Satellite data on such objects could provide us a key 
to the solution of this intriguing problem. 
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